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A NEW MERIT FUNCTION FOR EVALUATING THE FLAW 
TOLERANCE OF COMPOSITE LAMINATES 

by 

Martin M. Mikulas, Jr. and Rod Sumpter 


ABSTRACT 

Advanced composite materials are providing notable weight savings in 
numerous aerospace and non-aerospace applications by virtue of their 
uniquely high strength-to-density ratio and high stiffness-to-density-ratio. 
However, their full potential for high performance has yet to be achieved 
because of the low tolerance of the "brittle" material to impact, holes, flaws, 
and localized discontinuities due to relatively low toughness. In the present 
paper a new and simple merit function is proposed for maximizing the 
strength performance of flawed composite laminates for different design 
requirements. This new merit function takes into account the flaw tolerance 
of a composite laminate as well as its stiffness. After development of the 
new merit function, an example is presented of its application to a practical 
laminate family. 


INTRODUCTION 

Advanced composite materials are providing notable weight savings in 
numerous aerospace and non-aerospace applications by virtue of their high 
strength-to-density ratio and high stiffness-to-density-ratio. However, their 
full potential for high performance has yet to be achieved because of the low 
tolerance of the "brittle" material to impact, holes, flaws, and localized 
discontinuities. Numerous studies have been conducted over the past fifteen 
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years to develop new failure theories and to experimentally characterize the 
flaw performance of composite materials. An important observation made 
in these studies was that flawed graphite/epoxy laminates tend to fail at 
strain levels around 0.0041. in other words, even though an unflawed 
laminate may be able to sustain strain levels greater than 0.01, the same 
laminate with even slight flaws could fail at strain levels on the order of 
0.004. Because of this observation, composite structures are presently 
designed such that their ultimate strain levels do not exceed this lower level 
associated with a flaw. 

The primary mechanisms contributing to this low ultimate strain level are 
low resin toughness and high local stresses due to stress concentrations. 
Over the past ten years, the materials community has been exploring ways of 
improving resin toughness and this work is leading to improvements in 
flawed-laminate ultimate strain capacity. The second area, associated with 
high stress concentrations, also offers a potential for improving the flaw 
tolerance of composite materials because of the large number of design 
variables available in the laminate. However, a major obstacle to exploring 
laminates with lower stress concentrations is the lack of an appropriate merit 
function for consistent evaluation of laminate performance. 

The design of a laminate for a specific application typically involves a 
delicate balance between material strength and structural stiffness. For 
example, a high-aspect-ratio aircraft wing will require specified levels of 
axial and torsional stiffness while maintaining a high level of strength, 
which may be distinct from that required for a low-aspect-ratio wing. On the 
other hand, the need for torsional stiffness may not be as high for a launch 
vehicle. For most orbital space structures the primary need is axial stiffness 
with lower requirements for torsional stiffness and strength. Thus, 
depending upon the application, the laminate requirements will be different. 

In the present paper a new and simple merit function is proposed for 
assessing the strength performance of flawed composite laminates for 
different design requirements. This new merit function takes into account 
the flaw tolerance of a composite laminate as well as its stiffness. A 
demonstration of the new merit function is presented for the flaw case of a 
large circular hole. Since this is only one condition that may be critical for a 
strength-designed laminate, other potential failure modes must be 
considered. For example, the critical design requirement could be associated 
with a crack, a small hole, a local flaw, or impact damage. In order to 
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determine whether improved laminates can be obtained for these conditions, 
the designs would have to be evaluated with an approach similar to that used 
for the large hole in the present paper. For these other failure modes, either 
a point stress failure criterion or a fracture mechanics approach would have 
to be used for evaluation. After development of the new merit function, 
examples are presented of its application to graphite/epoxy laminates. 

MERIT FUNCTION DERIVATION 

Composite laminates are unique in that they are made up of constituent 
lamina that can provide a wide spectrum of effective material properties. 
Consider, for example, the in-plane properties of an orthotropic laminate 

family made up of a (O m ,±0 n ,9O p ) set of lamina. By changing the 

percentage of each of the three sets of lamina, and by changing the angle 0, 
the four independent global orthotropic elastic constants, E x , 
E y ,G xy , and v xy , can be changed over a large range. However, for specific 

applications there exist specifications on the relative stiffnesses that highly 
constrain the design freedom available. Another requirement of extreme 
importance in the design of a laminate is its capacity to carry load in the 
presence of a flaw. The term flaw herin is used to denote any deviation from 
a perfect plate such as a hole, a crack, or impact damage. 

It is well recognized that the flaw tolerance of a laminate can change 
considerably depending upon its effective in-plane orthotropic stiffness 
properties. For example, the stress concentration factor for a circular hole is 

quite low in a ±45° laminate (~2) but is quite high in a unidirectional 
laminate (~6). A corresponding fact for these laminates is that the 

extensional stiffness E x for the ±45° laminate is very low, while that for the 
unidirectional laminate is very high. Between these extreme cases there 
exist many laminates, each of which has different values of the stress 
concentration factor k and in-plane stiffnesses. A logical question that 
remains to be answered is, which is best from a flaw tolerance aspect? 
Although low extensional stiffness laminates tend to have lower stress 
concentration factors, they carry less load by virtue of their lower stiffness. 
Thus, there could be an optimum laminate which carries the most load in the 
presence of a flaw. In this section a merit function is developed to enable 
quantification of the flaw tolerance of the wide range of laminate lay-ups 
available. 
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Sketch a.- Flawed composite plate. 

A schematic of the flawed composite laminate to be studied is shown in 
Sketch a. For this laminate the weight is written as: 

W = pBtL (1.1) 

where p is the weight density of the material, and B, t, and L, are the width, 
thickness, and length of the laminate, respectively. In this derivation it is 
assumed that the design value of the in-plane loading is equal to the value of 
loading that the flawed laminate can withstand before failure, which is 
written as: 


Nx(Design) — N x (pi aw ) (1*2) 

Further, it is assumed that the failure load of the flawed laminate is the 
ultimate failure load of the unflawed laminate multiplied by some flaw 
tolerance factor \j/. This relation is written as: 

^x(Design) = ^x(Flaw) = V^xCUltimate) (1*3) 

As will be discussed later the flaw tolerance factor \j/ is a reduction factor 
£ 1 that is assumed to account for various flaw effects such as holes, cracks, 
or impact damage. 
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The constitutive relations for plane-stress in an orthotropic medium are: 

1 
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In the present paper only a uniaxial loading is considered so that Ny = Nxy = 
0. Thus, from equation (1.4) the ultimate strength of the uniaxially loaded 
balanced and symmetric laminate is written in terms of the ultimate strain of 
the laminate as: 

Nx(Ultimale) = ^x^x(Ultimate) (1*8) 


Substituting this into equation (1.3) yields: 

^x(Dcsign) — V^x^x(Ultimate) (1*9) 

Equation (1.9) can now be solved for the required laminate thickness t, 
which in turn can be substituted into the laminate weight equation (1.1) to 
obtain: 




W = pBL 


\ 


vV E xy 


N 


x(Design) 


(Ultimate) 


(U0) 


A further assumption made in this derivation is that the ultimate failure 
strain, £ x (uitimate)> is not a function of laminate construction. This 
assumption should be relatively good for laminates with a large enough 
percentage of 0° fibers that the failure of the laminate is governed by the 
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failure strain of the 0° fibers (fiber-dominated failures). However, further 
exploration of this assumption is beyond the scope of the current paper. 
With this assumption, and since p, B, L, and N x (D es ign) are all constant for 

a given design, this equation shows that the minimum weight design for a 
flawed laminate occurs when the product \|/E X is a maximum. In the next 
section, a study is made of the variation of this newly proposed merit 
function \|/E X for one laminate family and one flaw condition. 

DESIGN OF A LAMINATE FOR IMPROVED FLAW TOLERANCE 
IN THE PRESENCE OF A LARGE CIRCULAR HOLE 

To demonstrate the application of the merit function \j/E x , the design of a 
laminate required to carry a specified load in the presence of a large circular 
hole is considered. The term "large" herein refers to hole diameters large 
enough that the laminate failure is governed by the stresses due to a stress 
concentration factor such as described in references 1 and 2. In references 1 
and 2, it was shown that for a finite-width graphite/epoxy laminate, failure 
under uniaxial loading can be closely predicted by using the stress 
concentration factor when the hole diameter is larger than five inches, 
independent of hole-diameter to plate-width ratio (see figure 5-b of 
Reference 2). For this case, the flaw tolerance factor is given by: 

v = i (2.1) 

k 

where k is the stress concentration factor for the laminate. Thus the merit 
function to be maximized is E x / k. In the present paper the merit function 
is demonstrated for the case of an infinite orthotropic sheet with a hole. 

Determination of In-Plane Elastic Constants and Stress Concentration 
Factor For an Orthotropic Sheet 

For a balanced laminate, the engineering in-plane orthotropic elastic 
constants can be determined using 



(3.2) 
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In equations (3.1) through (3.4) 



(3.3) 

(3.4) 


A ij = XQijtk (3*5) 

k=l 


and h = total laminate thickness. 

The Qy' s for the lamina te are the transformed lamina stiffnesses of classical 

laminate theory and are presented in Appendix A. For an axially loaded 
infinite orthotropic sheet with a circular hole (see reference 3) the stress 
concentration factor is written in terms of the engineering elastic constants 
as 


k = 1+ 




+ 



(3.6) 


Example Application of New Merit Function 

For a baseline laminate, the 48-ply graphite/epoxy laminate presented in 
reference 2 was chosen; For purposes of determining in-plane orthotropic 
properties only, the laminate can be considered to be composed of only three 
layers. These three layers each have thickness t^ and consist of all the 

0°plies, all the |±45°| plies, and all the 90° plies, respectively. The lamina 
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and laminate material properties for this 48 ply laminate are listed in Table 1 
where t* =t £ /h. 


Material Properties 

Lamina 

Laminate 

Longitudinal extensional modulus, (psi) 

19x106 

10.2x106 

Transverse extensional modulus, (psi) 

1.89x106 

5.05x106 

In-plane shear modulus, (psi) 

9.28x105 

2.93x106 

Major Poisson’s ratio 

0.38 

0.48 


Baseline laminate ply orientation: (±45 / 0 2 / ±45 / 0 2 / ±45 / 0 / 90) 2s 
(t ply = 0.0055", h = 0.264") 


Laminate layer 

Layer thickness ratio, t* 

to 

0.417 

t45 

0.50 

t% 

0.083 


Table 1: Elastic and Thickness Properties of Example Laminate 

By changing the relative values of the t/s and varying the angle of the 

±45° ply to some angle ±0, the flaw tolerance of the resulting laminates can 
be determined using the merit function Ex/k. All calculations were 
performed using equations (3.1) to (3.4), and equation (3.6) with the 
stiffness constants given in Appendix A. 

In this example, the dimensionless thickness t 90 was left constant at t 90 = 
.083, its value in the baseline laminate. The values of t 0 and t e were then 
varied from 0 to .917 while holding the sum t 0 + 1 0 = 0.917 a constant, and 
0 was varied from 20° to 60°. The effect that these changes have on the 
laminate stiffness Ex and stress concentration factor k is shown in Figures 1 
and 2 respectively. As expected, E x increases nearly linearly for increasing 
to- The stress concentration factor k, however, increases in a non-linear 
manner. The merit function, E x /k, is plotted for the same range of 
parameters in Figure 3. Also included on each figure is a horizontal line 
corresponding to the ordinate parameter for an isotropic laminate of the 
same constituent lamina. 
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Figure 3 indicates that maximum values of the merit function Ex/k do occur, 
and that there are values of the merit function greater than that for the 
baseline laminate. To study the characteristics of the laminates with high 
values of the merit function, four cases shown as dots and labeled 1 to 4 in 
the figures are investigated herein. The cases are described as follows 

1) For this case, the laminate angle was held constant at 45°, and the 
thickness t 0 was increased until the merit function Ex/k became a maximum 
at a value of t 0 = 0.7 (see Figure 3). The resulting increase in the merit 
function is about 10% for this case. However, as shown in Figure 4, the 
shearing stiffness decreased by 39%, which may not be acceptable for some 
applications. The corresponding Poisson's ratios for these laminates are 
given in Figure 5 for sake of completeness in studying the characteristics of 
die laminates. 

2) For this case the thickness t 0 was held constant and the angle 0 was 

reduced to 25°, resulting in an increase in the merit function of about 12%. 
Again the shearing stiffness decreased significantly, 28% in this case. 

3) & 4) For these cases both to and 0 were changed in such a fashion as to 
hold the shearing stiffness equal to that of the baseline laminate. For case 
(3), t 0 =0.35 and 0 = 35°, while for case (4), t 0 =0.24 and 0 = 30°. For both 
of these cases the increase in the merit function is about 5%. 

These increases in the merit function would represent corresponding 
decreases in the laminate weight for strength-designed laminates, as 
indicated by the weight equation (1.10). It should also be pointed out that 
the achievable increases in the merit function found in this paper are for 
thickness and angle variations away from nominal for the specific 48 ply 
laminate selected as an example. These percentage increases will change 
depending on the initial laminate. Thus, the performance of a flawed 
laminate using this new merit function should be maximized on a case-by- 
case basis. 

Isotropic Laminate Properties 

As a point of reference, the properties for an isotropic laminate are included 
on each figure. To accomplish this, the effective in-plane properties for an 
isotropic laminate were calculated using the graphite/epoxy lamina values 
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listed in Table 1. These resulting isotropic values are E = 7.73x106 psi, v = 
0.318, and k = 3. For this isotropic laminate, the value of the merit function 
is Ex/k = 2.6x106 psi. For comparison, the value of the merit function for 

the baseline case of the previous example was 3 . 1 x 106 p S i. Thus, the 
baseline orthotropic laminate has a flaw tolerance performance 20% better 
than an isotropic laminate. 

Merit Function Discussion and Implications 

Application of the newly proposed merit function \|/E X for flawed laminates 
was demonstrated using an infinite graphite/epoxy laminate with a circular 
hole. This laminate was a 48 ply (±45/0 2 /±45/0 2 /±45/0/90) 2s lay-up. 
Results are also given for a corresponding isotropic laminate. To provide 
insight into the characteristics of laminates which have improved flaw 
tolerance, parametric plots in figures 1-5 were made of E x , k, E x /k, G X y, 

and v X y for a general (o m ,±0 n ,9O p ) laminate. From the results of the 
parametric plots the following trends are observed: 

1) Although the stress concentration factor k increases with increasing 
laminate stiffness, the axial stiffness E x increases faster resulting in a higher 
value of the merit function Ex/k, (Figures 1, 2 and 3). 

2) Increases in the merit function tend to be accompanied by reductions in 
the laminate shearing stiffness G X y (Figure 4). 

3) Laminates with high values of the merit function also tend to have high 
values of the major Poisson's ratio v X y (Figure 5). 

The structural implication of these trends can be properly evaluated only in 
the context of integration of the laminates into a built-up structure for 
specific applications. 


CONCLUDING REMARKS 

In the present paper a new and simple merit function is proposed for 
maximizing the strength performance of flawed composite laminates for 
different design requirements. An example is also presented of using the 
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merit function to improve the flaw tolerance performance of a practical 
graphite/epoxy laminate. 

It is well recognized that the flaw tolerance of a laminate can change 
considerably, depending upon its effective in-plane orthotropic stiffness 

properties. For example, the stress concentration factor of a ±45° laminate 
is quite low (~2) while the stress concentration factor for a unidirectional 
laminate is quite high (~6). A corresponding fact for these laminates is that 

the extensional stiffness for the ±45° laminate is very low, while that for the 
unidirectional laminate is very high. Between these extreme cases there 
exist many laminates, each of which has different values of the stress 
concentration factor and in-plane stiffnesses. A logical question that 
remains to be answered is, which is best from a flaw tolerance aspect? 
Although low-stiffness laminates tend to have lower stress concentration 
factors, they carry less load by virtue of their lower stiffness. Thus, it seems 
that there could be an optimum laminate which carries the greatest load in 
the presence of a flaw. 

The present paper introduces the concept of a flaw tolerance factor \jf which 
accounts for the ultimate load reduction caused by a flaw. It is then shown 
that the minimum weight of a strength-designed laminate in the presence of 
a flaw will occur when the merit function, \j/Ex, is a maximum. Thus, the 
use of this merit function permits a rational selection of laminate lay-up for 
maximizing its flaw tolerance. 

An example of the use of this merit function is presented for the case where 
the laminate flaw is a large circular hole. This is the simplest flaw case to 
study because failure is well predicted by the laminate stress concentration 
factor. The results of the application of this merit function to other flaw 
types is beyond the scope of the present paper. 

To study the application of this merit function, a practical 48-ply 
(±45 / 0 2 / ±45 / 0 2 / ±45 / 0 / 90 )^ graphite/epoxy laminate was selected from a 
previous reference. The relative values of the thicknesses of the plies and 
the angle 0 were varied to obtain the maximum value of the merit function 
\yEx. The results from this study indicate that an increase of 10% could be 
obtained for the merit function with a corresponding decrease in laminate 
weight. It is also noted that this percentage is valid only for the particular 
laminate studied, and that the merit function would have to be applied on a 
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case-by-case basis depending upon the specific structural design 
requirements. 

As a further demonstration of the use of the merit function, its value was 
calculated for an isotropic laminate for comparison with the orthotropic 
(±45 / 0 2 / ±45 / 0 2 / ±45 / 0 / 90)^ laminate. Results of this study showed that 

the orthotropic laminate's flaw tolerance merit function was 20% greater 
than that of the isotropic laminate. Thus, the use of the merit function 
provides a simple, rational basis for discriminating among candidate 
laminates for a specific application. 
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APPENDIX A: Orthotropic Elastic Constant Expansion 

In this expansion, t* is the dimensionless thickness of each of the three 
layers of the laminate. For an orthotropic lamina the Qy* s are defined as: 


Qn = 


Hit 

1-Vi2 v 21 


Q22 = 


E 22 

I-V12V21 


V 12 E 2 2 

12 1-v 12 v 21 


Q66 = G 12 


The Qjj* s for the three-layer laminate considered in this paper are found by 
using the stiffness transformation equations from reference 4. 

For the 0° plies 

Q 11 =Qn 
Q 22 - Q 22 
Q 12 = Q 12 

066 = 066 

For the 90° plies 

Q 11 = Q22 

022= Qll 
Q 12 = Ql 2 
066=066 

For the ±0 plies 

Ql 1 = Ql 1 cos4 9 + 2 (Qi 2 + 2Qgg ) sin 2 0 cos 2 0 + Q 22 sin 4 0, 

Q 22 = Qi 1 sin 4 0 + 2 (Qj 2 + 2Q6g)sin 2 0cos 2 0 + Q 22 cos 4 0, 

0l2 = (Qll + Q 22 ~ 4Q 66 )sin 2 0cos 2 0 + Q 12 (sin 4 0 + cos 4 0 ) 

066 = (Ql 1 + 022 “ 2Qi2 “ 2Q66)sin 2 0cos 2 0 + Q66(sin 4 0 + cos 4 e) 
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Substituting these equations into equation (3.5) and simplifying gives the 
laminate stiffnesses, Aij, as: 


4u. = E >1 (to +?e c °s 4 e) 

h 1 — V 12 V 21 ' 

+ — ^22 — (toQ + 1 6 (2 Vi 2 sin 2 0cos 2 0 + sin 4 0)) 

1 — V19V21 ' ' '' 


V12V21 
+4Gi2t0 sin 2 0 cos 2 0 


tZL = E H (r 90 +i e sin 4 8) 

h l-v 12 v 2 i v ' 

+—^ 22 — ^0 + te (2 vj 2 sin 2 ©cos 2 0 + cos 4 ojj 

44G 12 t0 sin 2 0cos 2 0 


■^12. = — ^11 (to sin 2 0 cos 2 0) 

h l-v 12 v 2 i V ' 


+- — ^22 — |vi2(tQ + tyo) + 10 sin 2 0cos 2 0 + v^^sin 4 0 + cos 4 ©jj 


■V12V21 
—4Gi2t0 sin 2 0 cos 2 0 


a 66 E ll 

h l-v 12 v 2 i 

+ 


|t e sin 2 0 cos 2 ©j 


— ^ 22 — [t 0 s i n 2 9 cos 2 0 — 2v 12 t 0 sin 2 0cos 2 g) 
1-Vi 2 v 2 i v ’ 

+Gi2^o + tpo - 2t 0 sin 2 0cos 2 0 + t 0 (sin 4 0 + cos 4 0 


For an orthotropic lamina, the minor Poisson’s ratio v 2 j can be determined 
using the reciprocity relation (reference 5). 

V21 _ V12 

^22 ®11 
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k, Stress concentration factor 



Figure 2 - Stress concentration factor for example laminate 
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Figure 5 - Major Poisson's ratio for example laminate 
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